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Abstract. We have analyzed six OB supergiants and one giant cov ering spectral types from 03 to Bl in the Galactic OB as- 
sociation Cyg OB2 by means of an updated version of FASTWIND (Santolaya-Rey et al., 1997) that includes an approximate 



treatment of metal line blocking and blanketing. This large coverage in spectral type allows us to derive a new temperature 
scale for Galactic O supergiants that is lower than the one obtained by using pure H-He models, either plane-parallel and 
hydrostatic or spherical with mass-loss. The lower temperatures are thus a combined effect of line blanketing and the large 
mass-loss rates. In some cases, the newly derived effective temperature is reduced by up to 8 000 K. Changes are larger for 
earlier stars with large mass-loss rates. As a consequence, luminosities are modified as well, which results in a lower number of 
emerging ionizing photons and reduces the mass discrepancy. Although there are still significant differences between spectro- 
scopic and evolutionary masses, we do not find any obvious systematic pattern of those differences. We derive mass-loss rates 
an d the correspond ing wind momentum-luminosity relation for the analyzed stars. Although consistent with previous results 
by Puis et al., 1996 for Galactic stars, our relation is better defined due to a reduction of errors related to stellar distances and 
points to a possible separation between extreme Of stars (Of 1- , Of*) and stars with more moderate morphologies. However this 
finding is only tentative, as the statistics are still scarce. 
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1. Introduction 

New winds are blowing from hot massive stars. The de- 



cent years ( 


Hubeny & Lanz, 1995; Santolaya-Rey et al., 1997 


Hillier & Miller, 1998; 


Pauldrach et al., 2001 


) employing im- 



proved numerical techniques and detailed atomic models of- 
fers the opportunity to derive more realistic stellar parameters 
that can be used for the research of astrophysical problems. 
Thus very recently Martins et al. (2002) have presented a new 
temperature scale for massive O dwarfs that lowers consider- 
ably the previous one given by Vacca et al. (1996) as a re- 
sult of strong metal line blanketing. The effect should be even 
larger for supergiants, where strong mass-loss should add to 
metal line opacity to yield even smaller effective temperatures. 
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First indications of this behaviour were obtained for exam- 
ple by Crowther & Bohannan (1997), Herrero et al. (2000) or 
Fullerton et al. (2000). More recently, Crowther et al. (2002) 
have presented an analysis of 4 supergiants in the LMC and the 
SMC with similar trends and Bianchi & Garcia (2002) obtain 
comparable results for 06-07 stars in the Milky Way. 

The stellar masses of O stars are still a subject of debate. 
A large part of the uncertainty originates from the so-called 
mass-discrepancy (Herrero et al., 1992): masses derived from a 
spectroscopic analysis using hydrostatic, plane-parallel model 
atmospheres are systematically lower than those predicted by 
non-rotating evolutionary models. Unified model atmospheres 
are expected to reduce the discrepancy, especially if new evo- 
lutionary models account for rotationally induced mixing in 
stars ( |Herreroet al., 2000| ; |Meynet & Maeder, 2000| ). To this 
end, considerable effort has been made in the last decade to ac- 
count for the possibility of He enriched and CNO contaminated 
atmospheres in early type stars. Efforts to improve the derived 
abundances have lowered them, but without completely ruling 
out He enhancements ( Villamariz et al., 2002| ). 
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On the other hand, the presence of a wind momentum- 
luminosity relation (WLR) for hot, massive stars offers the 
unique opportunity to derive stellar distances directly from the 
analysis of the observed spectra with an accuracy of about 10% 
(Kudritzki et al., 1999). For this purpose, the coefficients in the 
relation 

logiMV^R 5 ) w —,logL + C = x log(L/L e ) + logD 
a 

have to be calibrated as function of metallicity and spectral 
type. For O-stars, this has been done by Puis et al. (1996) for 
Galactic and MC objects. However, their WLR shows a large 
scatter, especially for non-supergiants. It is reasonable to as- 
sume that part of this scatter results from uncertainties in the 
individual stellar distances, as their sample was mainly taken 
from Herrero et al. (1992), aiming at targets with large appar- 
ent brightness. 

Following this idea, we started a programme to calibrate the 
WLR using a more homogeneous sample of OB supergiants, all 
of them belonging to the same OB association. From this ap- 
proach we expected to reduce the scatter in the derived WLR. 
The OB association chosen was Cyg OB2, and the optical ob- 
servations and first analysis using plane-parallel, hydrostatic 
models were presented in Herrero et al. (1999). The required 
wind terminal velocities were obtained from HST STIS obser- 
vations and are discussed in Herrero et al. (2001). 

By means of spherical, mass-losing models we will derive 
here an improved WLR for those Cyg OB2 supergiants that 
were observed with the HST. (One of the stars, Cyg OB2 #4, 
is actually a luminosity class III object, but we will refer to 
our sample as Cyg OB2 supergiants. When required, we will 
emphasize the difference). 

This effort is being complemented in our group by other 
studies to calibrate the WLR in the Local Group (see McCarthy 
et al., 1997, Urbanejaet al., 2002, Bresolin et al., 2002 forM31 
and M33). 

The remainder of this paper is organized as follows. In 
Sect. U we describe the characteristics of the code we have 
used, and in Sect. || we present an analysis of the 09 V star 
10 Lac as a test case. The analyses of our targets is given in 
Sect. |], and we end with a discussion of our results in Sect. |] 
and the conclusions in Sect. |J 

2. The code 

The analyses presented here have been performed by means of 
the latest version of FASTWIND (an acronym for Fast Analysis 
of STellar atmospheres with WINDs), a code which was orig- 
inally described by Santolaya-Rey et al. (1997) and was re- 
cently updated to include an approximate treatment of metal 
line opacity effects, i.e., metal line blocking and blanketing. 

The code follows the philosophy of using suitable physical 
approximations allowing for a fast computational time, thus en- 
abling us to calculate a large number of models (which are nec- 
essary to analyze stellar atmospheres with winds), while being 
realistic enough to preserve their value as a tool for determining 
stellar parameters. 



It assumes a /3-velocity law in the wind and calculates a 
consistent photospheric structure; the temperature structure is 
approximated using "NLTE Hopf functions" as described in 
Santoloaya-Rey et al. (1997); the coupling between the ra- 
diation field and the rate equations is treated within an ALI 
scheme, using local ALOs following Puis (1991). The program 
allows for a solution of the rate equations using either Sobolev 
or comoving frame calculations, but all calculations presented 
here were done using comoving frame. The formal solution 
to calculate the emergent profiles utilizes a radial micro-grid 
to account for the different scales involved and includes Stark 
broadening, which is a prerequesite for the analysis of O stars 
by means of H and He lines. 

The approximate treatment of metal line block- 
ing/blanketing will be described in detail by Puis (2002), 
in the following we will give only a brief summary. The basic 
philosophy to calculate the required NLTE metal opacities 
(bound-bound, bound-free and free-free) follows the line of 
reasoning given by Abbott & Lucy (1985), Schmu tz (1991), 
Schaerer & Schmutz (1994) and Puis et al. ( |2000| ), however 
applying significant improvements. In particular and most im- 
portant for realistic results, we have reformulated the equations 
of approximate ionization equilibrium (e.g., Puis et al. 2000, 
Eq. 43) to account for the actual radiation field (as function 
of depth) at all ionization edges (including those from excited 
levels) and employ a consistent coupling of rate-equations and 
mean intensity, in a way similar to the ALI approach, to avoid 
any kind of Lambda Iteration. 

The underlying atomic data base has been described by 
Pauldrach et al. (1998). In order to save computational effort, 
the resulting metal line opacities are averaged in a suitable way 
(mean of inverse opacties, in analogy to Rosseland means) over 
a frequency interval of order wind terminal velocity before the 
radiation transport is performed. Finally, flux conservation (and 
thus line blanketing) is obtained by adapting the NLTE-Hopf 
parameters in a consistent way. The method has been care- 
fully tested by comparing the results with up-to-date meth- 
ods, in particular with calculations performed with TLUSTY 
( Hubeny & Lanz, 1995 ) for the ca se of dwarfs (see also the 
next section) and with WMBasic ( Pauldrach et al., 2001 ) for 
the case of stars with winds. 

3. A test case: 10 Lac 

Before applying the modified code to the analysis of our Cyg 
OB2 sample, we have analyzed the 09V star 10 Lac as a test 
case. This star is well suited for calibration purposes because 
it is a luminosity class V star with a low projected rotational 
velocity, has already been used as standard by Herrero et al. 
(1992) and was later considered by Hubeny et al. (1998) to 
study the effects of line blanketing in plane-parallel, hydro- 
static atmospheres. 

Fig. [j] gives the fit diagram for the H/He spectrum of 10 
Lac. It is very similar to the corresponding fit diagram in 
Herrero et al. (1992), but now it is centered at a lower effec- 
tive temperature (and a slightly lower gravity). A comparison 
of both diagrams also reveals that the dispersion around the fi- 
nal model is now smaller. Therefore, the present error box has a 
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width of only 1 000 K instead of 2 000 K, as the one in Herrero 
et al. (1992). It is reassuring that all lines, including Hell A 
4200 that could not be fitted by Herrero et al. (1992), lie within 
the error box. 

The stellar parameters of 10 Lac we have derived here are 
T cfr = 35 500±500 K, log 3 = 3.95±0.10 and e= N{ ^\ Rc) = 
0.09±0.03, N(X) being the abundance of element X by number. 
The only significant difference compared to the results from 
Herrero et al. (1992) is the effective temperature, now lower by 
2 000 K. This is in complete agreement with the result obtained 
by Hubeny et al. (1998). These authors estimated a temperature 
of 35 000 K for 10 Lac using TLUSTY, a plane-parallel, hydro- 
static, line-blanketed model. 

Fit diagrams have the drawback that they only give the best 
possible fit for the chosen constraints. They rely on interpola- 
tions and sometimes (when using EWs) do not account for the 
profile shape. The actual final fit may still be poor. Fig. || shows 
the line spectra for our final model of 10 Lac. Good agreement 
is found for all lines, although a few details are not perfectly 
reproduced. In particular, the core of He II A 4200 is too weak, 
which is also true for the forbidden component of He I A 447 1 . 
Besides this, however, the final fit is perfectly consistent in all 
other aspects. 

Our result also agrees with the study of Martins et al. 
(2002), who recentl y found that using pure H/He models (as 
Herrero et al., 1992 did) results in an effective temperature 



using their EWs. Our results from CMFGEN completely agree 
with those from FASTWIND. 

The reason for the lower temperatures derived is twofold 
(see also Martins et al., 2002). On the one hand, the radiation 



scale for O dwarfs hotter by 1 500 to 4 000 K compared to using 
line-blanketed models. 

These authors have derived a ne w effective temperatur e 
scale for O dwarfs using CMFGEN ( ^illier & Miller, 1998| ), 
a spherical code including mass-loss and blanketing. In their 
T e ff scale, 09V stars are located at 33 000 K. However, their 
scale is calibrated using the equivalent width (EW) - spectral 
classification relations of Conti & Alschuler ( 1 97 1 ) and Mathy s 
(1998), and thus we have to compare our result with Martins et 
al. using the classification of 10 Lac in this system. 

Conti & Alschuler (1971) have classified 10 Lac as 08 III 
and not as a luminosity class V star (which have been con- 
sidered by Martins et al., 2002 ), although the star lies just at 
the border between both luminosity classes. The luminosity 
class III is mainly due to the low EW quoted by Conti & 
Alschuler ( 1 97 1 ) for the He I A4 1 43 line . Other EWs quoted by 
these authors are consistent with our observations which show 
a much larger value for this line, resulting in a luminosity class 
V within their classification scheme. Thus, we conclude that 
10 Lac should be classified as 08 V in the system of Conti 
& Alschuler (1971) (whatever the reason for the low EW in 



Conti & Alschuler, 1971 was) 



The effective temperature in the Martins et al. scale for 08 
dwarfs lies between 36 000 and 35 000 K, in perfect agreement 
with our result. Therefore we regard our result as fully consis- 
tent with recent findings using more elaborate but also more 
expensive fully blanketed NLTE model atmospheres. 

Nevertheless we have analyzed 10 Lac with our version 
of CMFGEN that includes a consistent photospheric structure 



field which is backscattered due to the additional opacity pro- 
duces a larger (E)UV radiation field. On the other hand, due to 
line-blanketing the electron temperature rises in photospheric 
regions. Both effects favour a higher ionization degree at lower 
effective temperatures, compared to unblanketed models. This 
effect can be clearly seen in Fig. |3|. 

The analysis of 10 Lac gives us an indication of what we 
can expect when introducing metal line opacity (namely lower 
effective tempe ratures), either in the form of traditional line- 
blanketing (as |Hubeny et al., 1998) , |Martins et al., 2002| and 



(Najarro, 2002) and have fitted the line profiles instead of only 



references therein) or including the line-blanketing via adapted 
Hopf-parameters (as here). Although with our method we do 
not force flux conservations, in all models calculated here the 
flux is conserved to better than 3% at all depths, where the re- 
maining deviations have no impact on the emergent fluxes and 
profiles. 

Our analysis of 10 Lac also gives us an idea of the error 
bar we can expect for the stellar parameters. For a resolution 
of 8 000, a S/N of 200 and a projected rotational velocity of 40 
km s -1 , the estimated errors are ±500 K in T c g, ±0.1 dex in 
logg and ±0.03 in e. For what follows we shall note here that 
this analysis does not give us information about the mass loss 
rate or the [3 exponent in the wind velocity law, as the wind 
of 10 Lac is negligible. We obtain an upper limit of 10~ 8 M 
yr _1 , but the fit has the same quality for any value below that 
limit. (The fit presented here was performed with M= 10 -10 
M Q yr- 1 ). 

4. Analysis of the Cyg OB2 supergiants 

We make use of o ptical and UV spect r a that have been pre - 
sented elsewhere ( ^errero et al., 1999 , Herrero et al., 2001 ). 
The optical spectra have been newly rectified, resulting in a 
less pronounced bump in the 4630-4700 A region, which, how- 
ever, does not affect our present analysis. Basic data, together 
with previously determined parameters adopted here are listed 
in Table [[], while the new parameters determined in this work 
are provided in Table || Note that the gravities in that table have 
been corrected for the effect of cen trifugal acceleration (in the 
same way as in Herrero et al., 1992 ). This effect is small (mod- 
ifying only the last digit of the entries in our table) and the 
gravities actually used in our calculations are always given by 
the nearest lower value ending with "0" or "5" in the second 
decimal (i.e., a corrected value of 3.52 corresponds to a model 
value of 3.50). The corrected gravities have been used for the 
calculation of stellar masses. 

Errors adopted for the parameters given in Table ^ have 
been taken from the errors obtained for Cyg OB2 #7 (see be- 
low) or from those for 10 Lac, depending on the stellar param- 
eters, the spectrum quality and the fit conditions. A summary 
of all errors is given in Table § 

For Cyg OB2 #7, errors have been estimated "by eye" from 
a microgrid of models around the final one. At fixed (3 (the ex- 
ponent of the velocity law), this resulted in ±1000 K in T e g, 
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Fig. 1. Fit diagram to the H/He spectrum of 10 Lac using FASTWIND plus approximated line blocking/blanketing. The filled 
circle marks the location of the chosen model, and the box around gives the size of the adopted error box 



t° 10 dex in log g, ±0.05 dex in e and to°l dex in M. (3 is de- 
termined from the form of the H Q wings, and its uncertainty is 
estimated to be ±0. 10. This has an influence on the derived ef- 
fective temperature and mass-loss rates. Therefore, the adopted 
error for T e g has been increased to ±1500 K and that for M to 

+0.10 
-0.15- 

The errors for radii and masses depend on the error in the 
absolute magnitude. This is assumed to be ±0. 1 from the work 
by Massey & Thompson (1991). (The influence of the error in 
T e g on the stellar radius is marginal and has been neglected 
here). For the analysis we have adopted a microturbulence of 
10 km s _1 . Tests indicate that this parameter is of no relevance 
for the results presented here. 

In the following we will comment on the individual anal- 
yses. Further discussions about the individual stars were pre- 
sented in Herrero et al. (2001). 



Cyg OB2 #7 The final fit to the observed H/He spectrum of 
Cyg OB2 #7 is shown in Fig. Q The good fit to the H Q profile 
is accompanied by a much poorer fit to the other two Balmer 
lines (and also to H$, not displayed here). This behaviour re- 
produces the one found in Herrero et al. (2000): for stars with 



strong winds we could not obtain a consistent fit for all Balmer 
lines at a given mass-loss rate. The situation has improved with 
the new version of our code, but the discrepancy still reaches 
a 30% effect, by which the mass-loss rate has to be reduced 
(from 10~ 5 to 7.7 xl0~ 6 ) in order to fit Hp and H 7 . The 
other stellar parameters are not affected by this modification, as 
the fit to the other lines does not change. We adopt the mass- 
loss rate indicated by H Q as this line is much more sensitive 
to changes in this parameter and there is good general agree- 
ment between the mass-loss rates from H Q and radio fluxes 



(Scuderi et al., 1998). Our result also supports this conclusion, 



as the mass-loss rate derived here agrees with the upper limit of 
1.5xl0~ 5 M© yr" 1 quoted by Bieging et al. (1989) (1.6 x 10~ 5 
if we use our values for distance and Voo) in case the star is a 
thermal emitter (the authors classify the object among the prob- 
able thermal emitters). However, in addition to the error quoted, 
M could be lowered by an additional 20% if we would adopt 
the value indicated by H 7 . 

The singlet He I lines (AA 4387, 4922) give a poor fit to 
the observed spectrum, partly due to the difficulties in the con- 
tinuum rectification. Therefore we have given a low weight to 
these lines when determining the stellar parameters. However, 
since these lines react strongly to changes in stellar parameter, 
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Fig. 2. The fit to the H/He spectrum of 10 Lac using FASTWIND with approximated line blocking/blanketing. The ordinate gives 
the relative flux values. Note that the scales are different for each line. See text for comments. 



it is always possible to find a reasonable fit within the error 
box. The error in e is larger than for other stars and e itself 
is not well constrained towards higher He abundances, because 
the already large He abundance produces a saturation effect. 

The derived T c g is much lower than the one obtained by 
Herrero et al. (2000) using the same code as here but without 



line-blocking. The derived luminosity is also lower by more 
than 0.2 dex, as the radius does not change very much. The 
reason can be seen in Fig. |[ where we compare the emergent 
energy of two models for Cyg OB2 #7. The first model (the 
dash-dotted line in the figure) is the one adopted here, with a 
T e g of 45 500 K, a radius of 14.6 Rq and metal line opac- 
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Fig. 3. The He ionization fractions in two models with T c g= 35 500 K, log<?= 3.95 and e= 0.09, with approximated line- 
blocking/blanketing (solid line) and without (dashed line). At typical photospheric line formation depths the ionization degree 
of the model including metal opacities is larger. Thus, lower temperatures are required in these models to reproduce the He 
spectrum. 

Table 1. Cyg OB2 supergiants studied in this work. All numerical identifications are taken from Schulte (1958). Magnitudes 
have been adopted from Massey & Thompson (1991), as well as spectral types, except for Cyg OB2 #11 and =#4, that are taken 
from Walborn (1973). S/N values have been estimated from the rms of the continuum at several wavelength intervals. Velocities 
are given in km s ~ 1 . 



Went 


V 


Spectral 


M v 


S/N 


Vi sini 


Voo 




mag. 


Type 










7 


10.55 


03If 


-5.91 


140 


105 


3080 


11 


10.03 


05 If+ 


-6.51 


190 


120 


2300 


8C 


10.19 


05If 


-5.61 


195 


145 


2650 


8A 


9.06 


05.5 1(f) 


-7.09 


135 


95 


2650 


4 


10.23 


07 m«f)) 


-5.44 


230 


125 


2550 


10 


9.88 


09.51 


-6.86 


145 


85 


1650 


2 


10.64 


Bl I 


-4.64 


195 


50 


1250 



ity. The second model (solid line) corresponds to the model 
adopted by Herrero et al. (2001), with a T cS of 50 000 K, a ra- 
dius of 14.8 i?0 and pure H/He opacities. Both models give the 
same optical luminosity and thus fit equally well the observed 
visual magnitude of Cyg OB2 #7. We additionally show in the 
figure the CMFGEN luminosities calculated with the same pa- 
rameters and conditions (dashed lines). The good agreement 
supports our approximated treatment of the metal line opacity. 

The reason for the similarity in derived radii is that the 
emergent flux is strongly blocked in the UV by the metal line 
opacity and thus emerges at higher wavelengths, including the 
optical. Therefore, at lower temperatures we obtain larger op- 
tical fluxes for models that include metal line opacity than for 
pure H/He models of the same temperature. 

The radius needed to fit the observed visual magnitude is 
then significantly smaller for models with metal line opacities 
(compared to unblanketed models at the same T e ff)> however 
roughly similar to the "old" value derived from unblanketed 



models at higher T e g. In consequence, the reduction in lumi- 
nosity is mostly due to the change in the effective temperature. 
Note, however, the dramatic difference of the ionizing fluxes in 
the (E)UV. 

The helium abundance derived here is even larger than the 
one obtained by Herrero et al. (2000), although the error bars 
overlap significantly. Cyg OB2 #7 is thus confirmed as the only 
star in our sample for which we derive an enhanced He abun- 
dance. The (3 value we obtain (0.9) is slightly larger than the 
one obtained from the UV analysis presented in Herrero et al. 
(2001), a behaviour which has been found already in previous 
investigations (e.g., Puis et al., 1996[ ). However, our mass-loss 
rate is similar to the one derived in that work, resulting again 
from the fact that the optical fluxes are similar. 

Finally, the log g value obtained here is slightly larger than 
the one obtained by Herrero et al. (2000). This results in a 
spectroscopic mass of 39.7 M , to be compared with an evo- 
lutionary mass of 67.4 M©. Accounting for maximum errors, 
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Table 2. Results obtained using FASTWIND plus line blocking/blanketing. Effective temperatures have been derived from the 
He ionization equilibrium. Gravities include the correction for centrifugal acceleration, e is the He abundance by number relative 
to H plus He; R, L and M are given in solar units. M s , M ev and M are the spectroscopic, present evolutionary and initial 
evolutionary masses, the two latter from the models by Schaller et al. (1992). M is given in solar masses per year, and logMWM 
means the logarithm of the Modified Wind Momentum rate, Mv^ R 5 
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29.0 


3.11 
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30.3 


1.00 


3.09e-6 


5.77 


43.1 


43.7 


48. 


29.248 


2 


Bl I 


28.0 


3.21 


0.09 


11.3 


1.00 


6.92e-8 


4.85 


7.5 


18.1 


19. 


27.263 



Table 3. Errors adopted for the parameters given in Table ||, in the same units as in that table. If only one value is provided, it 
should be preceeded by the ± sign. See text for comments. 
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Type 
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logs 
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logR 
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logM 


logL 


M s 
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Mo 
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7 
11 

8C 


03 If* 
05If+ 
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1.5 
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+U.15 
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0.02 
0.02 


0.10 
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+U.1U 
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+0.10 
-0.15 
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0.10 
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+ 17 
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+ii 

-8 
+ 10 
-9 
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+ii 
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+U.12 
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8A 


05.5 1(f) 
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0.10 
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29 


+4 
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-6 
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-6 
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2 


Bl I 
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0.10 


0.14 


3 


3 


3 


0.13 



masses of 56.7 and 59.4 M , respectively, are possible, still 
not overlapping. However, we have to remember the large He 
abundance derived for Cyg OB2 #7. This is an indication that 
this star might be evolved or be affected by rotational mixing. 
Taken together, there is no (clear) evidence that evolutionary 
and spectroscopic masses really disagree. 



Walborn et al. (2000) and Walborn et al. (2002) have re 



cently studied Cyg OB2 #7 and HD 93 129 A. They conclude, 
from a comparison of their spectra, that Cyg OB2 #7 has 
to be cooler than HD 93 129 A, and in fact HD 93 129 A has 
been reclassified as the prototype of the new 02 If* class. 
Taresch et al. (1997) have analyzed this latter star and obtained 
an effective temperature of 52 000±1000 K, based on the 
N v/N IV/N III ionization equilibrium. This is similar to the 
Herrero et al. (2000) value and again much higher than the 
temperature we obtain here for Cyg OB2 #7. Both stars display 
simultaneously N V, N IV and N III lines in their spectra, and 
thus we would not expect a very large temperature difference. 
Clearly, a cross-calibration of He and N blanketed temperature 
scales for the earliest stars is an urgent task. 

Cyg OB2 #1 1 The fit to Cyg OB2#l 1 is given in Fig. g It 
shows the same problems as the fit to Cyg OB2 #7, except for 
the fit to the He I singlet lines, which are again affected by con- 
tinuum rectification problems. However their depths relative to 
the depressed local continuum are now well predicted. 

Therefore we adopt the same errors, except for e, for which 
we adopt ±0.03. It is interesting that in spite of the extreme 
Of character of both stars and the similarity of the problems 



found, we do not derive an enhanced He abundance for Cyg 
OB2 #11. This star also shows the same trend as Cyg OB2 #7 
towards cooler temperatures and lower luminosities, but now 
the spectroscopic and evolutionary masses (73.0 and 58.1 M , 
respectively) invert their usual ratio. When considering the for- 
mal errors presented in Table ||, the mass ranges of both stars 
overlap significantly. 

The f3 value we have used is again 0.9. The mass loss rate 
we derive is consistent with the upper limit given by Bieging et 
al. (1989) (1.4xl0~ 5 M yr" 1 , or 1.2xl0~ 5 using again our 
values for distance and wind terminal velocity). 



Cyg OB2 #8C The fit to Cyg OB2#8C is presented in Fig. ^ 
The only problem is a serious failure in the predicted Hell 
A4686 line (which is not used in the fit procedure). To fit this 
line one had to increase the mass-loss rate by at least a factor 
of 1.7, although we note that the observed line lies at top of 
a broad emission feature that we cannot reproduce. The same 
comments as for Cyg OB2 #1 1 apply for the He I singlet lines. 

is not well constrained from the wings of H Q , as these are 
in absorption. We have adopted the same value as for Cyg OB2 
#7 and #1 1 (0.9), as well as the same errors. 

The resulting mass-loss rate is consistent with the upper 
limit quoted by Bieging et al. (1989) of 8.8 x 10~ 6 M Q yr -1 . 
The gravity is large for a supergiant (a lower gravity is pro- 
hibited by the Balmer line wings), but we find good agreement 
between the spectroscopic and evolutionary masses. 
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Fig. 4. As Fig. 0, however for CygOB2 #7. See text for comments 



Cyg OB2 #8A The final fit can be seen in Fig. |. The fit to 
He II A4686 is problematic, but much less than for Cyg OB2 
#8C, while the same comments apply for the He I singlet lines. 
Errors are the same as for Cyg OB2 #8C. The temperature is 
cooler and the luminosity lower than quoted in Herrero et al. 
(2001). However, the spectroscopic mass is again very large 



(90.5 M ), larger than the evolutionary one (78.4 M ), but 
with significant overlap when considering the errors. 

The mass-loss rate we derive here is nearly a factor of two 
lower than the one given in Herrero et al. (2001). Note, that the 
latter was not derived from spectrum analysis, however was 
calculated from the luminosity (believed to be larger at that 
time) and the Galactic WLR. 
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Fig. 5. Emergent energy of two models for Cyg OB2 #7, each one calculated both with FASTWIND and CMFGEN. The first pair 
of models corresponds to the one adopted here including metal line opacity and is represented by the lower dashed (CMFGEN) 
and dash-dotted (FASTWIND) lines. The second pair corresponds to the model adopted by Herrero et al. (2001) that did not 
include metal line opacity and is represented by the upper dashed (CMFGEN) and solid (FASTWIND) lines. We see that all 
models give the same optical flux, but very different UV fluxes (see text). We also see the good agreement between FASTWIND 
and CMFGEN, although individual strong UV resonance lines are not visible in the former due to the approximate method 
applied. 



Our new mass-loss rate in Table g (1.35±gJ|x 10 -5 M 
yr _1 ) agrees well with the radio mass-loss rate given by 
Waldron et al. (1998)1 (1.97xl0~ 5 ), and lies between the ex- 



treme values one would derive from the fluxes given by Bieging 
et al. (1989) (1.1-6.1 x 10~ 5 ) assuming free-free emission. 
Although Cyg OB2 #8A is a known non-thermal emitter, with 
variable radio flux and spectral index (Waldron et al., 1998; 



Bieging et al., 1989), our H Q mass-loss rate is of the same or- 



der of magnitude as the radio mass-loss rates and consistent 
with their lower limit. This con sistency contradicts the sug- 
gestion by Waldron et al. (1998) that the X-ray emission might 
originate from an X-ray source deeply embedded in the stellar 
wind, i.e., a base corona model scenario, which would imply a 
much lower mass-loss rate («1.5x 10 -6 M Q yr _1 ). 



loss rates available to compare with (which is an indication of 
a rather low value). The derived effective temperature is still 
cooler than in Herrero et al. (2001), although the differences 
begin to decrease. The evolutionary and spectroscopic mass 
ranges agree within the large error bars. 



Cyg OB2 #1 The fit to Cyg OB2 #10 is given in Fig. [L0| As 
for Cyg OB2 #4, the main difficulties appear in the fit of the 
He II A4686 and the He I lines, where from the three He I lines 
one is slightly too strong, the second slightly too weak and the 
third one fits well. The errors adopted are the same as for Cyg 
OB2 #4. The derived effective temperature is still lower than 
in Herrero et al. (2001), but the difference is only 2 000 K. The 
spectroscopic and evolutionary masses agree well. 



Cyg OB2 #4 The fit to Cyg OB2 #4 is presented in Fig. §. 
The predicted He II A4686 and the singlet He I line at A4922 are 
too strong in the core, although we note the large scale in the 
corresponding plots. The fit of He I A4387 is acceptable taking 
the normalization into account. j3 is again not well constrained 
from the H Q wings and we adopt a similar value as for the 
cooler stars in our sample ((3= 1). However, the influence of (3 
on the other stellar parameters begins to decrease and therefore 
we adopt the same errors as for 10 Lac. 

The mass loss rate is not well constrained towards lower 
values, because the profiles react only slightly. In this case, as 
also for the next two Cyg OB2 stars, there are no radio mass- 



Cyg OB2 #2 The final fit to the star is shown in Fig. |TT} The 
adopted errors are the same as for Cyg OB2 #8C (because of 
the lower S/N compared to 10 Lac), although again the mass 
loss rate is not well constrained towards lower values. Here, 
the derived effective temperature is hotter than in Herrero et al. 
(2001), but also the derived He abundance has decreased sig- 
nificantly. The masses, however, do not agree, with the spectro- 
scopic mass being lower than the evolutionary one by a factor 
of two. 

Herrero et al. (2001) indicated some problems with the stel- 
lar classification as Bl I and its Cyg OB2 membership, be- 
cause this would result in a rather faint absolute magnitude 
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Fig. 6. As Fig. 0, however for CygOB2 #1 1. See text for comments 



for its spectral class. However, the large reddening quoted by 
Massey & Thompson (1991) indicates that the star is proba- 
bly related to or lies beyond Cyg OB2. On the other hand, this 
large reddening is comparatively low when compared to other 
Cyg OB2 members, which additionally weakens the above ar- 
gument. (Cyg OB2 #2 has the fourth lowest reddening in Table 
7 of Massey & Thompson, 1991, who list a total 64 Cyg OB2 



stars. The stars with the three lowest values lie in the same re- 
gion of the association as Cyg OB2 #2). Thus, we adopt the 
absolute magnitude derived from the canonical distance to Cyg 
OB2 and assume a larger error, ±0.2 instead of ±0.1, which 
also doubles the error in the (logarithmic) radius. 
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Fig. 7. As Fig. EL however for CygOB2 #8C. See text for comments 



5. Discussion 



5. 1. The effective temperature scale for massive O 
supergiants 



The results obtained in the preceding section allow us to ad- 
dress four important aspects of the properties of massive OB With the inclusion of mass-loss, sphericity and line- 
supergiants: the effective temperature scale, the ionizing fluxes, blocking/blanketing we derived lower temperatures than those 
the mass discrepancy and the wind momentum-luminosity re- quoted in Herrero et al. (2001), where we have used plane- 
lationship. parallel, hydrostatic models without line-blocking for all stars 
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Fig. 8. As Fig. EL however for CygOB2 #8 A. See text for comments 



except Cyg OB2 #7, for which a spherical model with mass- 
loss, but without line-blocking was used. As we have covered 
various spectral types from 03 to Bl in our analysis, we can 
obtain a temperature scale for O supergiants. We note however, 
that Cyg OB2 #4 is actually a luminosity class III star, while 
#2 is probably of class II. 



For our earlier supergiants (03-07) we find temperatures 
that are 4000 to 8 000 K cooler than in Herrero et al. (2001), 
while for the 09.5 I supergiant we obtain 2000 K less and for 
the B 1 star the temperature is now hotter, but with a lower He 
abundance. These findings are in qualitative agreement with 
theoretical expectations (e.g., Schaerer & Schmutz, 1994) and 



with the recent temperature scale for O dwarfs presented by 
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Fig. 9. As Fig. 0, however for CygOB2 #4. See text for comments. 



Martins et al. (2002). These authors obtain lower temperatures 
for O dwarfs by 1 500^-000 K when including line-blanketing 
and sphericity. 

Our temperature scale is shown in Fig. ^2j together with 
the one from Vacca et al. (1996) for O supergiants. The ef- 
fective temperature scale from these authors has been mainly 
derived from analyses using pure H-He, plane-parallel, hydro- 



static models, like those performed by Herrero et al. (1992). It 
is thus not surprising that our new scale is cooler. 

In spite of the still low number statistics, we already appre- 
ciate some interesting features in Fig. [l2|. First, from 03 to 09 
we see a smooth temperature decline into which Cyg OB2 #4 
fits perfectly, despite its luminosity class III. From this smooth- 
ness is excluded Cyg OB2 #11, much cooler than the other 05 
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Fig. 10. As Fig. |[ however for CygOB2 #10. See text for comments. 



supergiant, Cyg OB2 #8C, and even cooler than the 05.5 su- 
pergiant, Cyg OB2 #8A. The main difference in their proper- 
ties is the extreme Of character of Cyg OB2 #11, which thus 
appears to be related to cooler temperatures. This is indicating 
to us that all spectral signatures have a significance in terms 
of stellar parameters, and thus a temperature scale using only 
spectral subtypes of O supergiants, with the various nuances in 



their classification scheme, will neccesarily be of limited accu- 
racy. 

We also see that both temperature scales converge towards 
later spectral types, until the Bl star, Cyg OB2 #2. This object 
has a temperature that does not seem to fit into the general be- 
haviour, although data are still too scarce to know whether this 
has any significance. 
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Fig. 11. As Fig. |[ however for CygOB2 #2. See text for comments. 



The new temperature scale and the lower luminosities will 5.2. Ionizing fluxes from O SUpergiantS 

have an influence on other aspects, e.g., on the emergent ioniz- , ,. . . 

r, , , , i • i i j j i Our treatment of U V metal line opacity is approximate (in the 

ing fluxes that are now much lower than in the older models, as , . r XTT ™, , 

, i . . c „ „„, sense that we use an approximate NLTE approach and suit- 

we have seen in the discussion of Cyg OB2 #7. rr rr 

ably averaged lines opacities), and we do not pretend to give 

a detailed description of the UV radiation field. However, the 

emergent fluxes, should be correct in an average sense (i.e., 

neglecting distinct spectral features, see again Fig. |^), in par- 
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Fig. 12. Our temperature scale for O supergiants (circles) com- 
pared with the scale by Vacca et al. (1996) (triangles). New 
temperatures are much lower, except for the relatively cool Cyg 
OB2 #2 (B II). Note that the entry at 07 is actually a luminosity 
class III star. 



Fig. 13. Left: ionizing fluxes (photons cm~ 2 s _1 ) at the stellar 
surface versus effective temperature. Right: ionizing luminosi- 
ties (photons s _1 ) versus spectral type. Dots correspond to the 
stellar analyses presented here and lines represent the calibra- 
tions by Vacca et al. (1996) for O supergiants. 



ticular concerning frequency integrated quantities. (Note that 
the differences bluewards from Hell 228 A are majorly due 
to different temperature structures in the outer wind.) Thus, we 
should obtain a rather correct description of quantities extend- 
ing over broad spectral regions, like the number of photons ca- 
pable of ionizing hydrogen. 

This number, of course, is of extreme relevance to stud- 
ies of Hll regions surrounding the stars. Vacca et al. (1992) 
have calculated ionizing fluxes from plane-parallel, hydro- 
static, LTE, line-blanketed Kurucz-models (Kurucz, 1992) and 
concluded, from a co mparison with more elabor ate models 
available at that time (Schaerer& Schmutz, 1994), that their 



ionizing fluxes should be reliable. That would indicate that 
line-blocking/blanketing is the major effect when calculating 
the UV continuum ionizing flux at a given effective tempera- 
ture. 

Our calculations seem to support this conclusion, but also 
indicate that for an application to H II regions further effects 
have to be considered. 

In Fig. [l3| a we have plotted our calculated H ionizing fluxes 
versus the effective temperature of the corresponding model. 
The dashed line gives the calibration by Vacca et al. (1996) 
for O supergiants (their Table 7). We see that the agreement 
is good. Even more, the entry departing slightly from the re- 
lation corresponds to Cyg OB2 #4, that has a luminosity class 
III. Note that the ionizing fluxes are basically independent of 
stellar radius. Thus, at given temperature the effects of metal 
line opacity remain the major ingredient which detemine the 
emergent UV flux, and the differences between NLTE and LTE 
models seem to be small regarding the integrated hydrogen 
Lyman fluxes. (This behaviour should become significantly dif- 
ferent when, considering, e.g., He II ionizing fluxes because of 
the extreme sensitivity of this ion to NLTE effects, in particular 
as function of mass-loss.) 

In Fig. pj| b we have plotted the H ionizing luminosities as 
function of spectral type. Here, the differences with the calibra- 
tion by Vacca et al. (1996) are apparent. They clearly originate 



from our new relation between spectral type and effective tem- 
perature (and our lower radii for most stars, compared to the 
calibrations by Vacca et al., 1996 ), i.e., the stellar luminosities 
at given subtype are now smaller. (This discrepancy between 
both plots, of course, results from the former inconsistency of 
calibrating spectral types via unblanketed models, however cal- 
culating the number of ionizing photons from blanketed ones). 
Furthermore, for one case (Cyg OB2 #8 A) we find that our ion- 
izing luminosities are larger as a consequence of the larger ra- 
dius in this particular case. Since the ionizing luminosity is the 
quantity which really matters for the ionization of H II regions, 
our results indicate that statistically there are fewer photons 
available compared to earlier findings, but also that individual 
cases have to be studied in detail, because they can depart from 
the general trend. 

5.3. The mass discrepancy 

The mass discrepancy is a term that was used by Herrero et 
al. (1992) to refer to the systematic difference between spec- 
troscopically determined stellar masses of OB supergiants and 
their evolutionary counterparts. The latter were always system- 
atically larger, well beyond the error bars. 

In Fig. [l4] we have plotted t he evolutionary mass (derived 
from the "classical" models by |Schaller et al. (199*2) without 
rotation) and the spectroscopic masses obtained in our analysis 
for the Cyg OB2 supergiants. We see that the situation is by 
no means satisfactory: despite the (very) large errors adopted 
for the spectroscopic masses, three of the seven stars still do 
not cross the one-to-one relation and for two other we find only 
a marginal agreement. In fact, only two out of the seven stars 
have masses that agree reasonably well. However, compared to 
previous diagrams the situation seems to have improved: there 
is no clear systematic trend any longer. Roughly the same num- 
ber of data points lie on each side of the one-to-one relation and 
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Fig. 14. A comparison of the spectroscopic and evolutionary 
masses for the Cyg OB2 supergiants. Although there are still 
serious differences, no obvious systematic trend is present. See 
text for details and discussion. 

the apparent scatter might be related to problems in the individ- 
ual analyses. 

Nevertheless it is somewhat too early to conclude that the 
atmosphere models with sphericity, mass-loss and metal opac- 
ity agree with the evolutionary models without rotation, and 
thus give the same answer in a statistical sense. First, the large 
scatter in Fig. |l4| still poses a question for the masses of both 
set of models; and second, the stars in Cyg OB2 have moderate 
projected rotational velocities and are very young, with only 
the earliest type exhibiting an enhanced He abundance. Thus, 
older or faster rotating stars may have masses that disagree even 
more when derived using different methods. 

This result indicates that we badly need both a calibration 
of present mass scales based on early type binaries (or any other 
reliable method), and CNO abundances for isolated massive 
stars. 



5.4. The Wind momentum-Luminosity Relationship 
in Cyg OB2 

The original purpose of our work was to obtain a better con- 
strained WLR for Galactic O stars by using objects belonging 
to the same association, and thus minimizing the scatter intro- 
duced by uncertainties in the relative distances. The distance to 
Cyg OB2 and thus the derived luminosities may be in error, but 
all determinations will be affected in a similar way. 

Fig. |l5] displays the WLR obtained for our Cyg OB2 sam- 
ple. Errors for Voo are taken from Herrero et al. (1991). The 
stars #8A, #8C, #4 and #10 follow a nice sequence with low 
scatter. Cyg OB2 #7 and #1 1 seem to lie above this sequence. 
This is interesting because these stars display the most extreme 
Of signatures in their spectra, which might be related to an 
ionization change in the wind that could result in a different 
line-driving force or to clumping effects that would produce 
an overestimation of the mass-loss rate. Cyg OB2 #2 seems 
to lie below that relation, which is consistent with the results 
by Kudritzki et al. (1991), who found a different WLR (with a 
lower offset) for the winds of early B supergiants, compared to 



the O-star case. Our observed relation also agrees well with the 
theoretical WLR derived by Vink et al. (2000) (the dashed line 
in Fig. 15). 

We have investigated whether a reduction in the outer min- 
imum electron temperature (that the model is allowed to reach, 
see Santolaya-Rey et al., 1997, Sect. 3.1) might result in lower 
mass-loss rates for Cyg OB2 #7 and #11. While we usually 
assume that T m iii= 0.75 T ff (a typical value for OB stars), cal- 
culations by Pauldrach (private communication) indicate that 
this value can reach values as low as 0.40 T c ft in extreme 
cases. This effect has usually no influence on an analyis as per- 
formed here, where almost all considered lines are formed in 
a region with temperatures beyond this minimum. "Only" H a 
(and He II A4686) might be affected in cases of extreme mass- 
loss, which is the reason that we have investigated here this 
question. Nevertheless, in all considered cases the resulting re- 
duction of the derived mass-loss rate was less than 20%. Thus, 
even accounting for this uncertainty we cannot bring the posi- 
tion of these two stars into agreement with the WLR defined by 
the other four O stars. 

However, taking into account the error bars, our data are 
still compatible with a unique relation including all seven stars. 
This is also shown in Fig. |l5[ where we display two dif- 
ferent regressions, one including only Cyg OB2 #8A, #8C, 
#4 and #10, and one including all seven stars, respectively. 
Interestingly, the relation obtained by including only the O- 
supergiants, not shown for clarity, is still marginally compati- 
ble with the position of the B-supergiant. 

In Fig. |l6| we show a comparison of our data points with 
those obtained by Puis et al. (1996), updated by Herrero et al. 
(2000) for some entries. We see that both sets compare well, 
the most important difference being the fact that the scatter in 
our data points is lower. However, at present it is not possible 
to conclude whether this is a real improvement (as we have 
expected) or an artefact of the low number statistics. We note 
that the first possibility would imply that stars with extreme Of 
characteristics follow a slightly different relation than normal 
O stars, although we could not determine the reason for the 
relatively high position of Cyg OB2 #7 and #11. 

In Table |] we compare the coefficients we have obtained 
for the WLR (by weighted least-squares fits) using different 
samples (see above) as well as those quoted by Kudritzki & 
Puis (2000) in their Table 2, and the coefficients of the theo- 
retical relation provided by Vink et al. (2000). A correct error 
treatment would require to take into account errors in both axes 
and their correlation. This treatment is not s imple and will be 
presented elsewhere ( Markova & Puis, 2002 ). For the scope of 
the present paper we have adopted as error an average of the 
errors obtained when considering only those in the ordinate 
values (an underestimation) and when considering the errors 
in abscissa and ordinate as uncorrected (an overestimation). 
Using both estimates for the errors, we have calculated the cor- 
responding regression, where the resulting values for slope and 
offset turned out to be only marginally different. The final val- 
ues quoted in Table ^| have been obtained from a straight aver- 
age of these results. 

We see that in spite of the visual agreement found in 
Fig. |l^, our results differ significantly from those given by 
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Fig. 15. The wind momentum-luminosity relationship ob- 
tained for the Cyg OB2 supergiants. The dotted line gives the 
regression obtained by including all stars, while the solid line 
gives the one obtained by including only Cyg OB2 #8A, #8C, 
#4 and #10. The dashed line very close to the solid one is the 
theoretical relation by Vink et al. (2000). The entry with the 
highest luminosity is Cyg OB2 #8A and those close together 
are Cyg OB2 #7 and #11, the stars with extreme Of character. 
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Fig. 16. A comparison of our results for Cyg OB2 OB- 
supergiants (filled dots) with the results by Puis et al. (1996) 
and Herrero et al. (2000). 



Kudritzki & Puis (2000), in particular concerning the slope 
of the relation. Kudritzki & Puis (2000) found a slope of 
1.51 ±0.1 8, while, when excluding the B supergiant, we obtain 
values from 1.74±0.24 to 1.92±0.22. In addition, our values 
for the vertical offset vary between 19.27 and 18.30, below the 
value of 20.69 of Kudritzki & Puis (2000). Although the error 
bars allow for marginal agreement, the conclusion is that our 
relation is steeper and thus a' eS (= 1/x) is smaller in our case. 
This indicates a different slope of th e line-strength distribu- 
tion function (see Puis et al., 2000 and [Kudritzki & Puis, 200C , 



Sect. 4.1), i.e., a larger number of weak lines. 

In contrast, we obtain very good agreement with the the- 
oretical relation by Vink et al. (2000). Their relation, both in 
slope and offset, lies right between our relation when con- 
sidering all supergiants or only the moderate O, Of stars, re- 



spectively. Our observations are thus compatible with a unique 
WLR for all O supergiants, but favour a separation of the ex- 
treme O, Of stars. 

6. Conclusions 

We have analyzed a homogeneous set of spectra of Cyg OB2 
supergiants with spectral types ranging from 03 to Bl. The 
analysis has been performed using a new version of the code 
presented by Santolaya-Rey et al. (1997) that includes an 
approximate treatement of metal line-blocking/blanketing. A 
test applied to the 09 V star 10 Lac resulted in an excel - 
l ent fit, at temperatu res lower than those obtained previously 



( |Hubeny et al., 1998 



( Herrero et al., 1992 ) and in agreement with other findings 



[vlai-tins et al., 2002| ). 
The fits to the supergiants are also satisfactory, except for 
the discrepancy between the fits of H a and other higher B aimer 
lines for stars with d ense winds (see discussion in Sect. 6.1 of 
frlerrero et al., 2000| ) and the pure absorption profile predicted 
for He II A 4686 in Cyg OB2 #8C instead of the P-Cygni profile 
observed. 

We obtain temperatures that are cooler than in our previous 
analyses, as a consequence of the effects of sphericity, mass- 
loss and metal-line blocking/blanketing. Thanks to our cov- 
erage in spectral type, we can construct a temperature scale 
for O-supergiants. This temperature scale is 4000 to 8 000 K 
cooler at early types than the one presented by Vacca et al. 
(1996), mainly based on pure H-He analyses. An important ad- 
ditional result is that at the same spectral subtype, Cyg OB2 
#8C and #1 1 have a difference of 3 000 K in their effective tem- 
peratures, reflecting their very different mass-loss rates. For the 
hotter stars, additional information based on a nitrogen temper- 
ature scale is needed, as the He lines become progressively ei- 
ther insensitive to changes in stellar parameters or simply too 
weak. 

Models including metal line opacities result in cooler tem- 
peratures and similar radii than those using pure H/He opac- 
ities, therefore giving smaller luminosities. Together with the 
blocking produced in the (E)UV, the corresponding ionizing 
luminosities may change drastically, which will have a large 
impact on studies of the surrounding regions. However, each 
particular case has to be treated individually. 

The results obtained indicate that Cyg OB2 stars are very 
massive. Only for two of the seven stars we found spectro- 
scopic or evolutionary masses below 40 M Q . Although there 
are still significant differences between spectroscopic and evo- 
lutionary masses, we do not find any obvious systematic pat- 
tern. Thus, at present we cannot conclude that there is a dis- 
crepancy between both sets of masses in Cyg OB2. However, 
this association is a very young one and the stars analyzed show 
only moderate projected rotational velocities: it might be pos- 
sible that rotationally induced effects simply have had no suf- 
ficient time to become apparent. It is remarkable that we found 
an enhanced He abundance for only one of the stars, Cyg OB2 
#7. 

Finally, we derived a new calibration of the wind 
momentum-luminosity relationship for O supergiants, includ- 
ing errors resulting from our analysis. Our data indicate that 
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Table 4. Coefficients for the wind momentum-luminosity relationship (see eq. 2) obtained in this work and taken from Kudritzki 
& Puis (2000) and Vink et al. (2000) 



logDo 



I ell 



Comments 



16.81±1.16 
18.30±1.82 
19.27±1.82 
20.69±1.04 
18.68±0.26 



2.18±0.21 
1.92±0.31 
1.74±0.32 
1.51±0.18 
1.826±0.044 



0.46±0.04 

52+ 10 

U.3»_o.09 

0.66±0.06 
0.548±0.013 



Using all Cyg OB2 stars 

Cyg OB2 #2 not included 

Inclu ding only Cyg OB2 #8A , #8C, #4, #10 

From kudritzki & Puis, 2000| (Table2) 

From |vinketal„ 2000| (Eq. 15) 



there might be a different relation for extreme Of and for mod- 
erate O, Of stars, respectively, perhaps indicating different ion- 
ization conditions or clumping in the wind. Considering only 
the moderate O, Of stars, we obtain a very low scatter in the 
relation, but this might reflect only the low number statistics. 
In fact, our still limited sample is also consistent with a unique 
relation including all stars, even the B supergiant. Clearly, more 
data are needed to disentangle whether the large modified wind 
momentum rate we derive for the extreme Of stars is a real, 
physical effect or just the result of our poor statistics. 
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